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(57) Abstract 



The invention provides a method for channel estimation in mobile radio communications which adaptively compensates for channel 
distortion on a block-by-block basis. Tht discrete-time channel impulse response is initially estimated with a given length and then 
truncated by using a sliding window. A cost function associated with the window is measured as the length and position of the window 
is adjusted over the channel impulse response and the cost function is compared with a threshold. The invention provides means to and a 
method for adaptively adjusting the length of the window Lt and the corresponding number of states in the equalizer, 2 L *r-l if appropriate. 
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CHANNEL ESTIMATION USING A SLIDING WINDOW TECHNIQUE 

This invention relates to channel estimation methods and apparatus in 
mobile radio communications and in particular to a receiver which adaptively 
5 compensates for channel distortion on a block-by-block basis. 

In digital mobile radio communications, transmission channels suffer 
from severe distortion due to frequency selective fading. In addition, channel 
characteristics are normally time-varying due to the relative motion of fixed and mobile 
stations. Therefore, in order to allow for reliable transmission, the receiver must be able 

10 to estimate and compensate for channel distortion on a block-by-block basis. 

Various channel estimation and channel equalization methods have been 
proposed in literature and are commonly used in practical systems such as mobile 
cellular communication systems employing the European wireless digital cellular 
standard "GSM". In most cases the receiver performs channel equalization on the 

15 received signal using Maximum Likelihood (ML) or Maximum A Posteriori (MAP) 
probability data estimation, based on the knowledge of the Channel" Impulse Response 
(CIR). Most practical systems employ training sequences to enable the CIR to be 
estimated before the equalizer star-up. Fast time varying, fading channels require the 
changing channel response to be tracked and adjusted dynamically by the receiver for 

20 the duration of the received signal. Tracking of the CIR may be performed by means of 
decision directed algorithms, where tentative decisions from the equalizer are used to 
update the initial CIR estimate. Examples of receiver systems which perform channel 
estimation and channel equalization may be found in the following articles: "Bit 
Synchronisation and Timing Sensitivity in Adaptive Viterbi Equalizers for Narrowband 

25 TDMA Digital Mobile Radio Systems", A. Baier, G. Heinrich and U. Wellens, Proc. 
IEEE Vehicular Technology Conference, June 1988, pp 377-384 [Reference 1]; 
"Correlative and Iterative Channel Estimation in Adaptive Viterbi Equalizers for 
TDMA Mobile Radio", ITG-Fachbericht No. 107, VDE Verlag, April 1989, pp 363-368 
[Reference 2]; "Simulation and Hardware Implementation of a Viterbi Equalizer for the 

30 GSM TDMA Digital Mobile Radio System", A. Baier, G. Heinrich, P. Shoeffel and W. 
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Stahl, Proc. 3" 1 Nordic Seminar on Digital Land Mobile Radio Communications, 
September 1988 pp 13.7.1 - 13.7.5, [Reference 3]. 

An equalizer of a given complexity can only cope with a certain delay 
spread of signal paths. Therefore, the estimated CIR is usually truncated and this 
5 truncated CIR estimate is used in a reduced state equalizer, of the type described in 
"Mobile Radio Communications" by R. Steels, pub. By Rentech Press and IEEE Press, 
pp 560-575. The reduced state equalizer selects an appropriate segment of the estimated 
CIR by sliding a window of length L T over the whole estimated response, calculating 
the energy contained within the window at each window position and identifying the 
1 0 window position where the energy within the window is maximum. 

However, if the actual CIR is either longer or shorter than L T samples, 
the above approach results in a degradation of the equalizer performance. In fact, as it is 
intuitive, assuming a channel response shorter than required produces a residual ISI 
which will affect the data estimation. On the other hand, assuming an impulse response 
15 longer than necessary means making use of one or more taps of the estimates channel 
which are essentially noise. There is thus a requirement for a reduced-state equalizer 
which provides improved equalizer performance. 

According to a first aspect of the invention there is provided a method of 
estimating channel impulse response in a communication system comprising: 
20 initially estimating the discrete-time channel impulse response with a 

given discrete length; and 

truncating the channel impulse response by using a sliding window, and 
characterised in that a cost function associated with the window is measured as the 
length and position of the window is adjusted over the channel impulse response and the 
25 said cost function is compared with a threshold. 

The overall performance of the equalizer is improved by neglecting the 
portion of the estimated CIR which is likely to correspond to just estimation noise. 

The invention provides means to and a method for adaptively adjusting 
the length of the window I r and the corresponding number of states in the equalizer, 
30 2 ir - 1 if appropriate. This may be achieved by computing a cost function for different 
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lengths and positions of the window, for example by measuring the out-of-window 
power for each window, and comparing that cost with a suitable threshold. 

Once L samples CIR have been obtained by standard means, the position 
of the truncated CIR of different lengths L T is then obtained by performing several 
5 maximum energy searches based on sliding windows of length Lp. For each window 
length, the out-of-window energy of the estimated CIR is compared with a suitable 
threshold. The threshold may be fixed, or it may be obtained from a previous estimate 
of the SNR (based on the non-truncated channel estimate). The invention may also be 
implemented by comparing the energy of the individual CIR taps with a suitable 
10 threshold. Again, the threshold may be fixed, or it may be obtained from a previous 
estimate of the SNR (based on the non-truncated channel estimate). 

A detailed description of a practical digital radio receiver is described 
below, by way of example, and with reference to the following figures in which: 

Figure 1 shows in outline a typical GSM digital radio receiver, 
1 5 Figure 2 illustrates the GSM "normal" burst format; 

Figure 3 illustrates a sliding window; 

A typical implementation of a digital radio receiver is shown in Figure 1 
and comprises an equalizer 10, deinterleaver 15 and a channel decoder 16. The 
equalizer includes a channel estimator 1 1 and data estimator 14. The channel estimator 
20 includes an initial channel estimator 12 and a channel truncator 13. 

The discrete-time received signal as received by the equalizer can be 

written as 

r(k)=Zb(k-l)h(l) + n(k) (1) 

where b(k) e{-l,l} are the transmitted data symbols or the (known) training sequence 
25 symbols, h(l) /= 1, 2, L represents the taps of the Channel Impulse Response (CIR) 

and n(k) indicates white Gaussian noise with zero mean and variance cr 2 . 

The equalizer must first estimate (12) the CIR h{l) , before beginning the 

data estimation process (14). In some cases, e.g. in a GSM standard receiver, the initial 

CIR estimation is commonly performed by means of correlative channel sounding, see 
30 for example the above mentioned references 1 and 3. The taps of the CIR estimate are 
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obtained by correlating the received signal r(k) with N=16 bits b(k) out of the 26 bits of 
training sequence, shown in figure 2, and represented by: 

A(/) = (l/A0Z^('>(/+0 < 2 > 

where 

5 i = 0, 1, ... N-l; 

The initial CIR estimate can generally be performed by ML channel 
estimation, giving: 

A = [A(0),A(l),...,A(I-l)] 7 ' = (B , B)- , B r r (3) 

where 

10 r = [r(0),r(l),...,r(Ar-l)] T 

B = [b(0),b(l),...b(#-l)] r 

It can be seen that, due to the good autocorrelation properties of the GSM 
training sequence B*B s NI and equation 2 is the particular case of the more general 
15 ML technique (equation 3). Once L samples (h(l)) of the CIR are estimated, an 
appropriate segment of samples (£ r (/) for L T <L) is selected (in the channel truncator 
13) by a maximum energy search based on a sliding window at length L T , see reference 
1, and as shown in figure 3. This approach produces a bit synchronisation with a 
resolution equal to the sampling period T and results in a reduced complexity equalizer 

20 with a trellis of 2^"' states. 

Assuming that for the training sequence bits b(k) 

= 5>(0K* +0 = W*.o , (4) 
substituting (1) in (2) yields 

h(l) = h(l) + v(/), 1=0, 1 L-l. (5) 

25 where 

v(0 = T7Z*('>(' + 0 i=0,l,...N-l, (6) 

Suppose that Lj taps of h (I) are selected as a CIR estimate hj (/) , when 
the actual channel response is shorter than L T bit intervals. From (5), it is clear that if 
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h(l) =0 for some index /, the corresponding sample h r (l) is just noise, and its use can 

only degrade the equaliser performance. An improvement may be achieved by 
adaptively adjusting the parameter Lj on a block-by-block basis, and consequently- 
adjusting the number of states 2 ij ~ x of the equaliser, in order to match the actual 
5 length of the channel impulse response. 

A suitable cost function, such as the square error of the truncated channel estimate, 
may be used to estimate length of the actual CIR and is given by: 

^(/ 0 .M = Sl*( / )-M')l 2 (7) 

(=0 

where h r (l) = h(l) for / e W = [/ 0 ,/ 0 + L r -l], and h T {l) = 0 for / «s W . Taking into 
10 account (5), the above quantity can be written as: 

J(/ 0 ,I 7 .) = £KOI 2 +SIW)| 2 (8) 

Therefore, choosing the value of iloXf) which minimise the cost (7) implies a trade- 
off between reducing the estimation noise and representing the relevant ISI terms. 
From (5) and (8), and assuming jA(/K(/)«0, the following equation is finally 

lew 

15 obtained: 

J(/ 0 ,L,) = 2X|v(/)| 2 +X|/«/)| 2 (9) 

/elf W 

However, in comparing the cost associated with different window lengths Lj, the 
above approach requires an accurate estimate of the squared value of the estimation 
error v(k) as defined in (6). From this, as an alternative to (9), one can consider the 
20 following process. Assume that, for each window length Lj, the best window position 
Iq is determined by a maximum energy search. Once a value of lo is associated with 
each length Lj, the best Lj can be chosen by comparing the out-of-window energy 
Y) I 2 corresponding to the different values of Lj with a suitable threshold. In the 

case where the out-of-window taps of the estimated channel do not contain any 
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relevant ISI terms, using (5) and (6) and by denoting with E L and £ ir the energy of 



5 Therefore, for each burst, the quantity El-Elt is computed and the validity 

of (10) is checked for different values of Lt. The lowest Lr satisfying the 

inequality (10) is an estimate of the length of the actual channel response, 

and may be then used both as a length of the CIR and for setting the 

number of states in the equaliser. 
10 The effectiveness of the invention has been asses by computer simulation 

for the case of a GSM receiver. The GMSK transmitted data symbols are obtained 

from the source bits by a rate 1:2 convolutionally encoding and interleaving , according 

to the GSM specifications for the TCH/FS class lb bits. At the receiver, the channel 

estimation is performed by using N=16 bits out of the 26-bits training sequence 

15 midambleofthe GSM normal burst, according to (2). The original L = 8 taps estimate 
are then truncated to Lj< 6 samples according to the rule (9) or (10). For data 
estimation, we employ a symbol-by-symbol Max-Log-MAP equalizer on the ISI trellis 
with state complexity 2^ . The equalizer soft-output data are deinterleaved and then 
fed to a symboI-by-symbol-Max-Log-MAP convolution channel decoder. 

20 Two channel profiles from "Digital cellular telecommunications system 

(phase 2+): Radio transmission and reception." ETSI GSM 05.05 (version 5.2.0), July 
1996 are considered. The first is the GSM typical urban area profile (TU50) with a 
relative speed of 50 Km/h and a delay spread of about 5jisec. The second propagation 
condition is the GSM typically hilly terrain profile (HT100), with a relative speed of 

25 lOOKm/h and a delay spread of about 20jisec. Considering the GSM symbol interval T 
= 3.69 jxsec, and taking into account that the GSMK pulse duration is about 3T, the 
above propagation conditions correspond to the two opposite cases where the invention 
provides a significant (TU conditions) and a negligible (HT conditions) performance 
improvement respectively. 



h(l) and h T (l) respectively, for signal-to-noise ratios 



CL\KD\ l )l<r 2 >7 



one has 



El'Elt= (L'L t ) (a2/N) < (L-L T ) (E LT /yN). 



(10) 



WO 99/56440 PCT/GB99/01313 

-7- 

To measure the performance of the different channel estimation 
algorithms, the means square error between the true and estimated channel responses are 
considered. In the case of the TU50 conditions, both strategies (9) and (10 produce a 
gain at the equalizer output and at the channel decoder output. For the HT100 profile, 
5 no significant improvement is achieved and the BER is similar to a conventional 
receiver. 
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Claims 

1. A method of estimating channel impulse response in a 
communication system comprising: 

initially estimating the discrete-time channel impulse response with a 

5 given length; and 

truncating the channel impulse response by using a sliding window, and 
characterised in that a cost function associated with the window is measured as the 
length and position of the window is adjusted over the channel impulse response and the 
said cost function is compared with a threshold. 
10 2. A method as claimed in claim 1 characterised in that the cost function 

is computed for each different length and/or position of the window. 

3. A method as claimed in claim 1 characterised in that the out-of- 
window power is used as the cost function. 

4. A method as claimed in claim 1 characterised in that the length and 
15 position of the window given the best result, when compared with the threshold, is 

chosen as the channel impulse response estimate. 
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